I-V characteristics of the vortex state in MgB 2 thin films 
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The current-voltage {I-V) characteristics of various MgB2 films have been studied at different 
magnetic fields parallel to c-axis. At fields (j,oH between and 5 T, vortex liquid-glass transitions 
were found in the I-V isotherms. Consistently, the I-V curves measured at different temperatures 
show a scaling behavior in the framework of quasi-two-dimension (quasi-2D) vortex glass theory. 
However, at /j,oH > 5 T, a finite dissipation was observed down to the lowest temperature here, T = 
1.7 K, and the I-V isotherms did not scale in terms of any known scaling law, of any dimensionality. 
We suggest that this may be caused by a mixture of a band vortices and it band quasiparticles. 
Interestingly, the I-V curves at zero magnetic field can still be scaled according to the quasi-2D 
vortex glass formalism, indicating an equivalent effect of self-field due to persistent current and 
applied magnetic field. 



PACS numbers: 74.70.Ad, 74.25.Qt, 74.25.Sv 
I. INTRODUCTION 



Since the discovery of the two-gap superconductor 
MgB2 in 2001, 1 the mechanism of its superconductiv- 
ity and vortex dynamics has attracted considerable in- 
terests. The two three-dimension (3D) 7r bands and 
two quasi-two-dimension (quasi-2D) a bands in this sim- 
ple binary compound seem to play an important role 
in the superconductivity, 2 as well as the normal state 
properties^! 4 - The two sets of bands have different en- 
ergy gaps, i.e., about 7 meV for the a bands, and about 
2 meV for the tt bands. 5 ' 6 And the coherent length of 
the 7T bands is much larger than that of the a bands 2 . 
Many experiments have demonstrated that the 7r-band 
superconductivity is induced from the a-band and there 
is a rich evidence for both the interband and intraband 
scattering. Owing to the complicated nature of super- 
conductivity in this system, its vortex dynamics may ex- 
hibit some interesting or novel features. Among various 
experimental methods, measuring the current-voltage (7- 
V) characteristics at different temperatures and magnetic 
fields can provide important information for understand- 
ing the physics of the vortex state. Up to now, the 
transport properties of MgB2 have been studied on both 
polycrystalline bulk samples 7 and thin films 8 . In both 
cases, the I-V characteristics demonstrated good agree- 
ment with the 3D vortex glass (VG) theory. This was 
partially due to the limited magnetic fields in the experi- 
ment. In addition, it has been shown that the properties 
of MgB2 are very sensitive to the impurities and defects 
introduced in the process of sample preparation, and the 
vortex dynamics must be influenced, too. Therefore, it 
is necessary to investigate the vortex dynamics in high 



quality MgB2 epitaxial thin films and to reveal the in- 
trinsic properties of the vortex matter in this interesting 
multiband system. In this paper, we present the I-V 
characteristics of high-quality MgB2 thin films measured 
at various temperatures and magnetic fields. The vortex 
dynamics in this system is then investigated in detail. 



II. EXPERIMENT 

The high-quality MgB2 thin films studied in this work 
were prepared by the hybrid physical-chemical vapor de- 
position technique 9 on (0001) 6i7-SiC substrates. All the 
films had c-axis orientation with the thickness of about 
100 nm. Fig. Q] (a) shows the 9-26 scan of the MgB 2 film, 
and the sharp (000Z) peaks indicate the pure phase of the 
c-axis orientation of MgB2. In order to show the good 
crystallinity of the film, we present in Fig. |Tfb) the same 
data in a semilogarithmic scale which enlarges the data 
in the region of small magnitude. It is clear that, besides 
the background noise, we can only observe the diffraction 
peaks from MgB2 and the SiC substrate, i.e., there is no 
trace of the second phase in the film. The c-axis lattice 
constant calculated from the MgB2 peak positions was 
about 3.517 A(bulk valueA: 3.524 A). The (j) scan (az- 
imuthal scan) shown elsewhere^ indicated well the six- 
fold hexagonal symmetry of the MgB2 film matching the 
substrate. The full width at half maximum (FWHM) 
of the 0002 peak taken on the film in 9-20 scan [MgB2 
0002 peak in Fig. [2(a)] and oj scan [rocking curve, shown 
in Fig. [He)] is 0.15° and 0.39°, respectively. The scan- 
ning electron microscopy (SEM) image in Fig. [ljd) gave 
a rather smooth top surface view without any observable 
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FIG. 1: (a) X-ray diffraction pattern of the MgB2 film on a 
(000f) 6-ff-SiC substrate in the 9-29 scan, which shows only 
the 000? peaks of MgB2 in addition to substrate peaks, in- 
dicating a phase-pure c-axis-oriented MgB2 film, (b) The 
semilogarithmic plot of the 9-29 scan, (c) The rocking curve 
of the 0002 MgB 2 peak, which shows the FWHM of about 
0.39°. (d) The SEM image of the MgB 2 film, which shows 
the smooth surface without obvious granularity. 



grain boundaries, which suggested that the film had a 
homogeneous quality. Ion etching was used to pattern a 
four-lead bridge with the effective size of 380 x 20 /im 2 . 
The resistance measurements were made in an Oxford 
cryogenic system Maglab-Exa-12 with magnetic field up 
to 12 T. Magnetic field was applied along the c axis of 
the film for all the measurements. The temperature sta- 
bilization was better than 0.1% and the resolution of the 
voltmeter was about 10 nV. We have done all the mea- 
surements on several MgB2 films, and the experimental 
data and scaling behaviors are similar; so, in this paper, 
we present the data from one film. 

In Fig. [2l we present the resistive transitions (R-T re- 
lations) of a MgB2 thin film measured at various mag- 
netic fields in a semilogarithmic scale. The current den- 
sity in the measurement was about 500 A/cm 2 , much 
smaller than the critical value for low temperatures, 
10 6 A/cm 2 iS. It can be determined from Fig. [2] that 
the sample had a superconducting transition tempera- 
ture of T c = 40.05 K, with a transition width of about 
0.5 K. Its normal state resistivity was about 2.45 /iflcm 
and the residual resistance ratio [= p(300 K)/p(42 K)] 
was about 6.4. The I-V curves were measured at various 
temperatures for each field, and then we got the electric 
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FIG. 2: Temperature dependence of resistive transitions for 
fioH = 0, 1, 3, and 6 T, with the current density j = 
500 A/cm 2 . 



field (E) and the current density (j) according to the 
sample dimension. The current density was swept from 
5 to 10 5 A/cm 2 during the I-V measurements. 



III. THEORETICAL MODELS 

In the mixed state of high-T c superconductors with 
randomly distributed pointlike pinning centers, a second- 
order phase transition is predicted between VG state and 
vortex-liquid state.— The I-V curves at different tem- 
peratures near the VG transition temperature T g can be 
scaled onto two different branches^ by the scaling law 
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The scaling parameter z has the value of 4-7, and v w 1- 
2; D denotes the dimension of the system with the value 
3 for 3D and 2 for quasi-2D 13 ; /+ and /_ represent the 
functions for two sets of the branches above and below 
T g . Above T g , the linear resistivity is given by 
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At T g , the electric field versus the current density curve 
satisfies the relationship 



E(J)\t=t. ~ j 



(*+!)/(£>-!) 



(3) 



In 2D superconductors at /j.qH = T, a Berezinskii- 
Kosterlitz-Thouless (BKT) transition was found at a spe- 
cific temperature Tbkt-— At Tbkt, E oc j 3 , which is a 
sign of the BKT transition. A continuous change from 
the BKT transition at zero field to a quasi-2D VG transi- 
tion, and then to a true 2D VG transition with T g = 0K 
was found in TIBaCaCuO film, 15 which shows a field- 
induced crossover of criticalities. 
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A 2D VG transition may exist in a true 2D system 

with T g = K, i.e., there is no zero-resistance state at 

any finite temperatures. The E-j curves can be scaled 
byie 
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where To is a characteristic temperature, i>2D ~ 2, and 
p > 1, while g is a scaling function for all temperatures 
at a given magnetic field. The linear resistance is given 
by 



p lin cxexp[-(To/Tf] 



(5) 



This 2D scaling law can be achieved in the very thin 
films^ or in highly anisotropic systems at high magnetic 
fields.!^ 



IV. EXPERIMENTAL RESULTS AND 
DISCUSSIONS 

A. Quasi-two-dimension vortex-glass scaling in the 
low-field region (poH < 5 T) 

The E-j characteristics have been measured at various 
magnetic fields up to f2 T. In Fig. [3] we show the typical 
example at poH = 1 T for (a) E-j curves and (b) the cor- 
responding p-j curves in double-logarithmic scales. It is 
obvious that when the temperature goes below some par- 
ticular value (this is actually the vortex-glass transition 
temperature T g according to following discussions), the 
resistivity falls rapidly with decreasing current density 
and finally reaches the zero-resistance state which is the 
characteristic of the so-called VG state. At the tempera- 
tures above T g , the resistivity remains constant in small 
current limit. The current density of 500 A/cm 2 used in 
p-T measurement shown in Fig. [5] lies in this linear resis- 
tivity regime from about 10~ 3 to 1 pQ, cm. Consequently, 
these data sets provide the basic information on scaling 
if the data are describable by the VG theory. 

The inset in Fig. 2] shows the data of the pu n versus 
(T — T g ) and the fit to Eq. (2). The data are the same as 
those shown in Fig.[5]for poH = 1 T, and the attempt T g 
value is 31.4 K. In this double-logarithmic plot, the slope 
of the linear fitting gives just the exponent of v(z+2 — D), 
and the determined value is 8.08 ± 0.05. In order to have 
reasonable values for v and z, the dimension parameter D 
needs to be chosen as 2, i.e., the investigated system has 
the property of quasi-2D, which is similar to the situation 
found in BiSrCaCuOj 13 ' 20 This is further supported by 
the VG scaling of the data at 1 T. As shown in the main 
frame of Fig. 01 the scaling experimental E-j curves form 
two universal branches corresponding to the data above 
and below T g (31.4 K) with v = 1.32 and z = 6.12. 
At very large current density or a temperature near the 
onset of superconducting transition, the free flux flow 
regime dominates and, hence, the data do not scale. The 
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FIG. 3: (Color online) (a) E-j characteristics measured at 
fixed temperatures ranging from 30 to 36 K for poH = 1 T. 
The increments are 0.30 K in the range from 30.00 to 31.20 K, 
and 0.25 K in the range from 31.50 to 34.00 K respectively, 
and finally 35 K on the top. The dashed line shows the po- 
sition of T g , and the symbols denote the segments that scale 
well according to the quasi-2D VG theory. The thin solid 
lines denote also the measured data, however, located outside 
the scalable range, (b) p-j curves corresponding to the E-j 
data in (a). The thick solid line in (b) denotes the voltage 
resolution of 10 nV. 
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Quasi-2D VG Scaling for |j H=1 T 
D=2, 1=31.4 K, v=1.32, z=6.12 
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FIG. 4: (Color online) Quasi-2D VG scaling of the E-j curves 
measured at 1 T. The inset shows a double- logarithmic plot 
of the temperature dependence of the linear resistivity. The 
dashed line is a guide for the eyes. 
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Quasi-2D VG Scaling for H=3 T 
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FIG. 5: (Color online) (a) Scaling curves of the E-j data mea- 
sured in 3 T based on the quasi-2D VG scaling theory. The 
inset shows a log-log plot of the temperature dependence of 
the linear resistivity, (b) VG scaling with another form of 
scaling variable j/(T \T - T g |' y(D_1) ) . 



symbols in the figure denote the range of the data well 
described by the scaling law. 

The situation at fioH = 3 T is similar to that at 
{.iqH = 1 T. As shown in Fig. [UJa), the determined pa- 
rameters are T g = 15.4 K, v = 1.17, and z = 6.58. Inter- 
estingly, the previous work on MgB2 film 8 indicated that 
the 3D VG scaling theory (D = 3) is a better choice in 
describing the I-V characteristics in this system, though 
this experiment was done at magnetic fields lower than 
1 T. Moreover, the I-V curves were demonstrated to scale 
well by using the argument of j/(T\T -T g \ ). The same 
conclusions were also drawn on the polycrystalline MgB2 
samples 7 . In order to clarify this issue, we also analyzed 
our data using the form suggested in Ref. 8. As shown in 
Fig.EJb), such a scaling with j/(T \T - T g \ v(D ~ 1] ) as the 
scaling variable is worse than that with j / \T — T g \ u< " D . 
Most importantly, the dimension parameter D is still re- 
quired to be 2 instead of 3 as proposed in Refs. 7 and 
8. This confusion can be easily understood in terms of 
the two-band superconductivity of MgB2. As we know, 
there are two types of bands contributing to the super- 
conductivity of MgB2, namely, the 3D it bands and the 
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FIG. 6: (Color online) p-j data at temperatures 1.7 K and 
4 K to 20 K with 2 K-step, for /i H = 6 T. Temperature of 
the isotherms increases from bottom to top. 



2D a bands. Therefore, the structure of the vortex mat- 
ter must be affected by both of them. Although the su- 
perconductivity of 7r bands, induced possibly by that of 
a bands, is much weaker, it provides a large coherence 
length with 3D characteristics in the low-field region. 
Therefore, the vortices in this system may be quasi-2D 
like and, at the same time, they can possess large cores 
characterized by the coherence length of the 7r band su- 
perfluid. In this sense, the quasi-2D scaling should be 
more appropriate than the 3D one. However, when a 
higher disorder is induced in the system, especially in 
the boron sites, the interband scattering gets stronger 
and the anisotropy decreases, which may lead to a 3D 
vortex scaling. In this case, a more rigid vortex line can 
be observed, especially at low fields^ The good quasi-2D 
scaling at 1 and 3 T demonstrated here suggests that the 
phase transition from VG to the vortex liquid in MgB2 
resembles that in the high-T c superconductors. Together 
with the data shown below, we can safely conclude that 
a vortex glass state with zero linear resistivity can be 
achieved in the low field region due to the presence of the 
finite superfluid density from the 7r bands. Regarding the 
VG scaling^ a principal requirement is a proper deter- 
mination of T g , namely the temperature with a straight 
log-E-logj curve in the low dissipation part. The toler- 
ance for T g variation is very small (about ±0.3 K). With 
an inappropriately chosen T g , the scaling quality dra- 
matically deteriorates and, simultaneously, the values of 
v and z quickly deviate from those reported above and 
those proposed by theory. This validates our analysis 
here. 
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FIG. 7: (Color online) The scaling of the E-j isotherms with 
quasi-2D VG model for poH — 6 T. The inset shows the 
deviation of the pn n vs T — T g (T g = 0) relation from the 
linearity in the double logarithmic scale. 



FIG. 8: (Color online) Attempted scaling of the data with 
2D VG model [Eq. (4)] for p H = 6 T. The inset shows the 
nonlinearity of the relationship between pu n and temperature, 
the solid line shows the theoretical curve of true 2D VG theory 
[Eq. (5)]. 



B. Anomalous vortex properties in high field 
region 

As shown in Fig. [21 when the magnetic field reaches 
6 T, no zero-resistance state can be observed down to 
the lowest temperature, here 1.7 K. Consequently, no 
VG transition exists above 1.7 K at this field, as shown in 
Fig.[H The shape of the curve at T — 1.7 K suggests that 
the resistivity goes to a finite value as the current den- 
sity approaches zeroi^ 3 . As shown in Fig. the pn n versus 
T — T g seriously deviates from linearity for any possible 
T g value, indicating the inapplicability of Eq. (2) in the 
present case. Correspondingly, the quasi-2D scaling law 
fails here. A natural explanation is that, with increas- 
ing field, the 3D supercurrent from it bands is seriously 
suppressed^* 2 ^ and the quasi-2D vortex structure trans- 
forms into a 2D-like one dominated by the a band super- 
fluid. In Fig. [51 we show our attempt to apply 2D VG 
scaling on the data. Surprisingly, this attempt also failed, 
even though this model has been successfully applied to 
the layered superconductors with large anisotropy (or 2D 
property) such as Tl- and Bi-based high-T c thin films at 
high magnetic fields ; 18 ' 19 

The most reasonable explanation for this anomaly is 
that the supercurrent contribution from the tt bands is 
much easier to suppress by the magnetic field than that 
from the a bands, since the gap in the n bands is several 
times smaller than that in the a band. We suggest that 
at high magnetic fields (above 5 T), a different vortex 
matter state is formed, composed of quasi-particles from 
the 7T bands and vortices formed mainly by the residual 
superfluid from the a band. The 7r-band quasiparticles 
diminish the long range phase coherence of the supercon- 
ducting phase, which leads to a finite dissipation. Once 
the long range superconducting phase coherence is de- 
stroyed by the proliferation of a large amount of these 



7r-band quasiparticles, neither 3D nor quasi-2D VG scal- 
ing is applicable. Such a mixed state is obviously diffi- 
cult to be simply described by any known scaling theory. 
Recently, scanning tunneling microscopy studies showed 
that the quasiparticles of the 7r bands disperse over all of 
the superconductor, both within and outside the vortex 
cores 25 , which strongly supports our arguments. This is 
the basis for the explanation of the nonvanishing vortex 
dissipation at high magnetic fields in a zero temperature 
limit found recently on MgB2 thin films. 23 



C. Self-field effect at /i H = T 

For a 2D layered superconductor in zerofield, the above 
mentioned BKT transition may exist and be reflected 
in the I-V characteristics 15 . In the present MgB2 sam- 
ples, we have not found any evidence of this transition 
in low magnetic fields which would be consistent with 
the quasi-2D (instead of 2D) configuration of the vor- 
tex matter. Moreover, both the E-j curves and the p-j 
curves (as presented in Fig. [5]) are similar to the situation 
of poH = 1 T. Considering the narrow transition width 
at zero-field, we did the measurement carefully with an 
increment of 0.05 K. Obviously, there is no E(j) curve 
which satisfies the E cx j 3 dependence, as expected by the 
BKT theory. Since the current can induce self-generated 
vortices, it might be interesting to look at whether the 
quasi-2D VG model applies here. 

Similar in Sec. IV A, we present py ln versus (T — T g ) 
in a double logarithmic plot. From this graph, we de- 
termined the exponent in Eq. 2 (as shown by the inset 
of Fig. fTUTa)). A good quasi-2D scaling was obtained 
with parameters T g = 39.94 K, v = 1.12, and z = 6.61, 
as presented in Fig. [TUT a). Using the parameters deter- 
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FIG. 9: (Color online) (a) E-j data at various temperatures 
from 39.7 K to 40.5 K with an interval of 0.05 K for fi H = 
T, the symbols denote the region, where the data are scaled 
(from 39.70 K to 40.30 K). Temperature of the isotherms 
increases from bottom to top. The dashed line shows the 
position of T g and the symbols denote the segments, which 
scale well according to the quasi-2D VG theory. The thin solid 
lines are also the measured data lying outside the scalable 
range, (b) p-j curves corresponding to the E-j data in (a). 
The thick solid line in (b) denotes the voltage resolution of 
10 nV. 



mined here, one finds a self-consistency with the value 
of v{z + 2 — D), as determined in fitting the linear re- 
sistivity [Eq. (2)]. Both the temperature dependence of 
piin and the scaling curves at /jqH = T are similar 
to the situation at small field \iqH = 0.1 T [shown in 
Fig. fTCT b)] and hqH — 0.5 T (not shown in this pa- 
per), except for the slight differences of the scaling pa- 
rameters. The scaling parameters including the ones at 
HqH = 0.5 T are listed in Table Q] It was proven that 
current and magnetic field exhibit analogous effects in 
suppressing superconductivity and generating quasipar- 
ticles in conventional superconductors^ Similarly, the 
current-induced self-field may lead to a similar effect in 
the vortex state as an applied magnetic field. Nonethe- 
less, the good agreement of this simple scaling law with 
the zero-field data is interesting and worth studying in 
detail. Moreover, the values of v and z for zero field are 
very close to those for /iqH = 1 and 3 T, indicating a 
similar vortex dynamics in the whole low-field region. 
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FIG. 10: (Color online) (a) Quasi-2D VG scaling of the data 
measured at T. The inset indicates a good linearity of the 
temperature dependence of the linear resistivity, (b) Quasi- 
2D VG scaling of the data measured at 0.1 T. The inset indi- 
cates a good linearity of the temperature dependence of the 
linear resistivity. 



TABLE I: Quasi-2D VG scaling parameters at different fields. 



VoH (T) 


T g (K) 


;/ 


z 


0.0 


39.94 


1.12 


6.61 


0.1 


39.28 


1.30 


6.08 


0.5 


35.95 


1.37 


6.42 


1.0 


31.4 


1.32 


6.12 


3.0 


15.4 


1.17 


6.58 



V. SUMMARY 

We have measured I-V curves on high-quality MgB2 
films at various magnetic fields and temperatures. At 
magnetic fields below 5 T including the zero field, the 
curves scaled well according to the quasi-2D VG theory 
instead of the 3D model, in good agreement with the 
multiband superconductivity of MgB 2 contributed from 
the strong 2D a bands and weak 3D 7r bands. At the 
fields above 5 T, the curves did not scale according to 
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any known VG scaling laws, accompanied by the disap- 
pearance of a zero-resistance state. Based on our result 
combined with recent tunneling experiments, a different 
vortex state was suggested, namely, a state where the 
vortices composed of the superfluid from the a bands 
move through the space filled with numerous quasiparti- 
cles from 7r bands. 
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